QCD Sum Rule Study of the Masses of Light Tetraquark Scalar Mesons 
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We study the low-lying scalar mesons of light u, d, s flavors in the QCD sum rule. Having all 
possible combinations of tetraquark currents in the local form, QCD sum rule analysis has been 
carefully performed. We found that using the appropriate tetraquark currents, the masses of a, n, 
fo and ao mesons appear in the region of 0.6 - 1 GeV with the expected ordering. The results are 
compared with that of the conventional qq currents, where the masses are considerably larger. 
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The light scalar mesons have been subject to inten- 
sive discussions for many years The expected mem- 
bers are cr(600), rc(800), / (980) and a (980) of flavor 
SU(3) nonet. The existence of cr(600) which is denoted 
as /q(600) in the Particle Data Group has been confirmed 
also by a model independent theoretical analysis 0. 

Yet, their nature is not fully understood [H, 0|. Be- 
cause they have the same spin and parity as the vacuum, 
J p = + , they may couple to many different modes. In 
the conventional quark model, they are 3 Pq state of qq. 
Their masses are, however, expected to be larger than 
1 GeV due to the p-wave orbital excitation. Further- 
more, the mass ordering in a naive quark mass counting 
of m u ~ rrid < m s implies m a ~ m ao < m K < m/ . 
In chiral models, they are regarded as chiral partners of 
the Nambu-Goldstone bosons (ir, K, r\, rj) [5[. Due to the 
collective nature, their masses are expected to be lower 
than those of the quark model. In chiral perturbation 
theories they are also described as resonances of meson- 
meson scattering, whose quark content is dominated by 
{qq) 2 @]. Recent discussions on the scalar mesons are 
then largely motivated by its tetraquark components. 

Tetraquark structure of the scalar mesons was pro- 
posed long ago by Jaffe with an assumption of strong 
diquark correlations 0. Due to the strong attraction 
in the scalar diquark channel, their masses are expected 
to be around 0.6 - 1 GeV with the ordering of m a < 
m K < mfaqa, consistent with the experimental observa- 
tion [1, 2, [§] . The form of qq-qq was also proposed @ . 
Recent lattice study also showed an indication of the 
tetraquark for f [lfj|. 

If such tetraquarks survive, they may be added to 
members of exotic multiquark states. The subject of the 
multiquarks is important, providing a new opportunity 
to study colored dynamics which has not been reached by 
conventional hadrons ll| . The problem is also related to 
the origin of hadron mass as we will briefly discuss here. 

In this letter, we would like to report the results of a 
systematic study of the masses of the tetraquark scalar 
mesons in the QCD sum rule. In the QCD sum rule 
one extracts hadron properties from two-point correla- 



tion functions computed by the operator product expan- 
sion (OPE) of QCD. The non-perturbative effects are 
then incorporated by vacuum condensates of QCD oper- 
ators. By comparing the theoretical correlation functions 
to phenomenological ones, one can determine phy sical 
quantities such as masses and coupling constants |l2lll3j|. 



QCD sum rule analyses become subtle for hadrons 
containing more quarks such as tetraquarks and pen- 
taquarks. Due to high dimensionality of the interpolat- 
ing field, one needs to calculate many terms in OPE with 
high dimension. At the same time, it becomes difficult to 
find a good Borel window with keeping the convergence 
of OPE. Another point we would like to address here 
is the proper choice of hadronic currents (interpolating 
fields). When the OPE has to be in any way truncated 
up to certain terms, unless the current is suitably cho- 
sen, the resulting OPE and the sum rule may not work. 
In general, there are several independent currents for a 
given hadron state. The optimal current can be searched 
by making their linear combinations, as has been tested 
recently for the exotic tetraquark state [Tl |. 



Let us construct tetraquark currents of J PC = ++ , 
by establishing the number of independent currents. Fol- 
lowing the method in our previous work we adopt 
the diquark construction, where the diquark and antidi- 
quark have the same color, spin and orbital symmetries. 
Therefore, they must have the same flavor symmetry, 
which is either symmetric (Gf (qq) gj 6f (qq)) or antisym- 
metric (3f (qq) (g 3f (qq))- Then we assume the ideal mix- 
ing in which only isospin symmetry is respected and the 
currents are classified by the number of strange quarks. 
Hence, denoting light u, d quarks by q, a currents are 
constructed as qqqq, k currents by qsqq and fo and ao 
currents by qsqs. 



Using the antisymmetric combination for diquark fla- 
vor structure, we arrive at the following five independent 
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currents 

S 3 = (ulCj 5 d b )(u a j 5 Cd[ - u b -y 5 Cd%) , 

T% = (uZCa^d b )(u a a^CdJ + u b a^C%), (1) 
A% = (ulC^d b )(u a ^CdJ + u b ^CdJ), 
P% = (ulCd b )(u a Cd[ -u b CdJ), 

where the sum over repeated indices (fj,, v, ■ ■ ■ for Dirac, 
and a,b, ■ ■ ■ for color indices) is taken. Either plus or 
minus sign in the second parentheses ensures that the di- 
quarks form the antisymmetric combination in the flavor 
space. The currents S, V, T, A and P are constructed 
by scalar, vector, tensor, axial-vector, pseudoscalar di- 
quark and antidiquark fields, respectively. The subscripts 
3 and 6 denote the color states of the diquarks (antidi- 
quarks) which are combined into the color representa- 
tion 3 C and 6 C (3 C or 6 C ), respectively. The currents 
for other members are formed by the following replace- 
ments in flj). k : (ud)(ud) — > (ud)(ds), fo : (ud)(ud) — > 
(us)(us) + (u •*-> d) , aa : (ud)(ud) — > (us) (us) — (u <->• d). 
More details will be discussed in a separate publica- 
tion 0. 

Using the tetraquark current 77 which is one of the cur- 
rents of ([I]) or their linear combination, we have com- 
puted the correlation function 

IW? 2 ) = i J d^xe^iOlT^xWmO), (2) 

in the OPE up to dimension eight, keeping the current 
quark masses m u , and m s finite. As the primary re- 
quirement, the spectral densities pope must be positive 
definite. If truncation of OPE is not good, it happens 
that they become negative sometimes. Using the disper- 
sion relation, it is equivalently written as 

iiope(^ rdsJ^L, (3) 

Jo s-q 2 -le 

where pope = IhiTIope/t'"- This is then equated to the 
integral over the physical (phcnomenological) spectral 
density p p hcn(s). 

The phenomenological spectral density is parameter- 
ized as a sum of one pole and continuum contributions. 
Assuming that the continuum part is approximated by 
the one of OPE (duality) [lj], 

Pphc„(s) - fS(s - M 2 ) + 9(s - S0 )-Inin O P E (s) , (4) 

where M and / are the mass and coupling constant of 
the physical state under investigation. In order to extract 
physical quantities efficiently by suppressing the contin- 
uum contribution, the Borel transformation is performed. 
Finally we arrive at the sum rule equation 

f 2 exp(-M 2 /M 2 ) = I ° ds P oMs)e- s/Ml , (5) 
Jo 



which determines the mass and the coupling constant. 
The mass M is a function of the two parameters Sq and 
Mb- They must be chosen to satisfy (1) rapid conver- 
gence of OPE, (2) sufficient amount of pole contribution 
and (3) weak dependence on so and Mb- These are im- 
portant in order to draw reliable conclusions [lij]. 

The use of the (^-function in (j4]) for the scalar mesons 
might be subject to criticisms, firstly because the ob- 
served scalar mesons have wide decay width. The in- 
clusion of a finite width in the QCD sum rule instead 
of using the ^-function is one option to take care of the 
effect of the decay width. We have performed such anal- 
ysis in the form of Gaussian and verified that still it is 
possible to reproduce experimental values of masses |15| . 

Secondly, the tetraquark currents is expected to couple 
strongly to two meson states. We argue that such two 
meson contributions can be computed separately from 
the short distance method of OPE. By applying the soft- 
pion theorem [13], the coupling to two meson states can 
be expressed by a double commutator with the axial 
charge Q 5 , (0 tt^tt 6 ) - (0|[Q 5 , [Qb,»?]]|0>. Since one 
commutator yields the factor qq, we find (qq) 4 altogether 
in the two-point correlation function. The dimension of 
this term is as high as twelve, which is beyond the present 
study where we compute up to dimension eight. Simi- 
larly, as shown in Ref. [l8( which also uses the method 
of QCD sum rule, the coupling of s — > pp (s stands for a 
scalar tetraquark and p pseudoscalar meson) is of higher 
order as proportional to (qq) 2 , consistently implying that 
the decay width of the s — > pp is of order (qq) A . 

As we will see shortly, the fact that the present QCD 
sum rule with the OPE up to dimension eight will yield a 
stable solution indicates that there is a significant com- 
ponent in the scalar meson state which couples to the 
tetraquark current without going through two mesons. 
In this case, we expect that the narrow resonance ap- 
proximation is reliable. The large decay width will then 
be explained through the coupling to two meson states 
in the form of high-dimension terms of OPE. In fact, we 
have performed a QCD sum rule analysis by using a peak 
of finite width, and found that the result does not change 
much. 

We have performed the sum rule analysis using all cur- 
rents and their various linear combinations. We have 
found that the results for single currents are not reliable, 
except for the tensor current Tg , due to either violation 
of positivity or insufficient convergence of OPE. In fact, 
we have found good sum rule by a linear combination 
of A% and V£: t\\ = cos9A% + sindV^ , where the best 
choice of the mixing angle turns out to be cot = 1/ 
For k, fo and ao, we have also found that similar linear 
combinations give better sum rules. 

The calculation of OPE is tedious but straightforward. 
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The results up to dimension eight are 
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fast [l5(. As Mb is increased, the pole contributions de- 
crease, but the resulting tetraquark masses are stable as 
shown in Fig. [T] 
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FIG. 1: Masses of the a (short-dashed), k (solid), fo and ao 
(long-dashed) mesons calculated by the tetraquark currents 
as functions of the Borel mass Mb, with so (GeV 2 ) as shown 
in figures. 
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The OPE for ao takes the same form as for fo- For tr, 
terms containing u, d quark masses m q are small. For 
instance, the term of m q (qq) of dimension four is about 
ten times smaller than the other term of (g 2 GG). For k, 
ao and /o, the terms containing strangle quark mass are 
important but those containing u and d quark masses are 
negligibly small. 

We use the following values of condensates [2(| 
H3: (qq) = -(0.240 GeV) 3 , (ss) = -(0.8 ± 0.1) x 
(0.240 GeV) 3 ,(g 2 GG) = (0.48 ± 0.14) GeV 4 , m u = 
5.3 MeV, m d = 9.4 MeV, m s (l GeV) = 125 ± 20 MeV, 
(g s qcrGq) = -A/ 2 x (qq), Af 2 = (0.8 ± 0.2) GeV 2 . Now 
let us discuss the feasibility of our QCD sum rule. 

The Borel transformed correlation functions are writ- 
ten as power series of the Borel mass Mb- Since the 
Borel transformation suppresses the contributions from 
s > Mfj, smaller values are preferred to suppress the 
continuum contributions. However, for smaller Mb con- 
vergence of the OPE becomes worse. Therefore, we 
should find an optimal value of Mb- We have found 
that M B ~ 0.4 GeV for a, 0.5 GeV for re and 0.8 GeV 
for fo and ao, where the pole contributions reach around 
50% for all cases, while the convergence is still sufficiently 



We have also searched the region where the tetraquark 
mass varies significantly less than the change in y / so. We 
have found such regions 0.5 < so(GeV 2 ) < 1.5 for a, 
1 < so < 2 for k, and 1.5 < So < 2.5 for fo and ao- 
In Fig. [51 we show So dependence of the masses in these 
regions. As we see, the mass is stable in a rather wide 
region of s . The Borel mass dependence of Fig. [T] are 
shown for the minimum values of so- 
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FIG. 2: Masses of the a (short-dashed), k (solid), fo and ao 
(long-dashed) mesons calculated by the tetraquark currents 
as functions of threshold value so, with Mb (GeV) as shown 
in figures. 



After careful test of the sum rule for a wide range 
of parameter values of Mb and sq, we have found re- 
liable sum rules, with which we find the masses m a — 
(0.6 ± 0.1) GeV, m K = (0.8 ± 0.1) GeV, to /o , Q() = 
(1±0.1) GeV . It is interesting to observe that the masses 
appear roughly in the order of the number of strange 
quarks with roughly equal splitting. 
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It would be interesting to observe from Eqs. ([6]) - ([8]) 
that the mass of the a is dominated by the gluon con- 
densate, while other condensates with m s also play a sig- 
nificant role for other masses. In fact in the SU(3) limit, 
where all quark masses and condensates take the same 
values, the three equations become identical. In partic- 
ular, in the chiral limit where all quark masses vanish, 
the masses of the scalar mesons are dictated only by the 
gluon condensate. This property was also observed in a 
recent publication j22|. 

Now for comparison, we have also performed the QCD 
sum rule analysis using the qq current within the present 
framework, although such works have been done be- 
fore [H HI 0, HI]. We have computed the OPE up 
to dimension six in this case, and have verified the pre- 
vious results. Namely, the masses of the qq mesons are 
considerably heavier than the masses of the tetraquark 
mesons. 

In conclusion we have found that the QCD sum rule 
analysis with tetraquark currents implies the masses of 
scalar mesons in the region of 600 - 1000 MeV with the 
ordering, m a < m K < m,f 0taal while the conventional qq 
currents are considerably heavier. Our conclusion has 
become rather robust, after we have tested all possible 
independent tetraquark currents and with their linear 
combinations. 

Our observation supports a tetraquark structure for 
low-lying scalar mesons. Somewhat non-trivial is that 
a large part of the mass is due to the gluon condensate 
rather than chiral condensate. This observation is in- 
teresting in relation to the ques tion of the origin of the 
mass generation of hadrons [27|. To test the validity of 
the tetraquark structure, it is also important to study 
decay properties, which is often sensitive to the struc- 
ture of wave functions. Such a tetraquark structure will 
open an alternative path toward the understanding ex- 
otic multiquark dynamics which one does not experience 
in the conventional hadrons. 
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